A dual order parameter from 
fundamentally colour charged matter* 

M. HOPFER, M. MlTTER, B.-J. SCHAEFER, R. ALKOFER 

Institut fur Physik, Karl-Franzens-Universitat Graz 
Universitatsplatz 5, 8010 Graz, Austria 



The center phase transition of QCD and of fundamentally charged 
scalar QCD at non-vanishing temperature is investigated within a Dyson- 
Schwinger approach. The temperature dependence of the scalar/quark 
propagator is studied with generalized boundary conditions. A novel order 
parameter for the center phase transition is established which still exhibits 
a considerable dependence on the scalar/quark-gluon vertex. 
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1. Motivation 

During the last years, intensive theoretical and experimental efforts 
towards a deeper understanding of the confinement and chiral symmetry 
breaking phenomena have been made. In these proceedings we focus on the 
QCD phase transition at non-vanishing temperature and the related center 
symmetry breaking^] by investigating so-called dual observables [31 HI [5j El 
[8] . Dual observables can be related to the spectrum of the Dirac operator 
and therefore comprehend information on both, the chiral properties and 
the confinement of quarks. They have first been introduced in lattice calcu- 
lations [3] but are also accessible via functional methods [5l|6j[71|8]. Within a 
Dyson-Schwinger approach [9j we study the center phase transition in QCD 
and in fundamentally charged scalar QCD. For the solution of the Dyson- 
Schwinger equation (DSE) of the scalar/quark propagator quenched lattice 
input for the gluon propagator has been used [7] . The main outcome of this 
investigation is a novel order parameter for the center phase transition |10j 
which still exhibits a considerable dependence on the used model for the 
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scalar /quark-gluon vertex. This requires a more detailed vertex analysis in 
future studies [TT] . 

2. Center Phase Transition and Related Order Parameters 

The DSE for the quark propagator at finite temperature reada^] 

f A'ijf 

S- 1 (p) = Z 2 S 1 (p)-Z 2 C F g 2 TY^ j^rS(k)T»(k,p-, q )D^( q ) (1) 

with the momenta p = (p,u) p ). The generalized Matsubara modes ojjd)) = 
(27rT)(n p + 4>/2tt) depend on the general boundary angle <fi 6 [0, 27r)|j The 
usual fermionic boundary conditions are obtained for <j) = tt. For the prop- 
agator we use the conventional ansatz .S -1 ^) = ij^ujpC(p) + ipA(p) + B(p) 
and for the quark-gluon vertex the temperature dependent model [6] 

^ A C(k) + C(p) , dv _ j A{k) + A{p) 

2 



di <Z 2 f«/i)ln[g 2 /A 2 + l]V^ 



I d2 + g 2 g 2 + A 2 \ 4-7T 

where and p denote the in- and outgoing quark momenta and q the gluon 
momentum. The first part of the vertex ^ consists of a Slavnov- Taylor 
motivated ansatz, whereas the logarithmic tail together with the anomalous 
dimension 25 = —18/44 ensures a running coupling-like behaviour of the 
vertex in the UV regime. The remaining model parameters d\ and d 2 are 
purely phenomenological, where in the following we take the values proposed 
in [7j . The system can now be solved by using the gluon propagator as input 
from quenched lattice results [7]. 

From the quark propagator ([I]), various order parameters can be ex- 
tracted. The ^-dependent quark condensate is given by 

= Z 2 N C T Y, J *r D S(p, (3) 

Up{<t>) 



Here, for the quark-gluon vertex the renormalization factor Z\f = Z2/Z3 has been 
used, while in Landau gauge the ghost-gluon vertex renormalization constant Z\ = \ 
[12] , Furthermore, the ghost renormalization constant Z3 is omitted as it cancels 
by a corresponding factor in the quark-gluon vertex model. The remaining constant 
Z2 is eliminated within a MOM scheme. Gp is the fundamental Casimir invariant, 
C F = 4/3 for N c = 3. 

A similar approach has been employed in [8], where QCD with ^-dependent boundary 
conditions for the quark fields has been referred to as QCD^. This corresponds to 
QCD at imaginary chemical potential. 



LightCone2012 printed on December 19, 2012 



3 



where the case <f> = ir corresponds to the ordinary quark condensate. With a 
Fourier transform of the previous condensate ^ the dual quark condensate 
[3] is obtained with respect to the winding number n 

"2tt ji 



En = J 2^ € ■ (4) 

In the following we fix n = 1. Si transforms like the conventional Polyakov 
loop pQ under center transformation and is therefore a suitable order pa- 
rameter for the center transition 0J El El E] , cf. also [S] . Alternatively, the 
scalar quark dressing function B, evaluated at vanishing external momenta, 
is also sensitive to center symmetry breaking [6j{^] However, if only the low- 
est Matsubara mode is taken into account as suggested in [6j no satisfactory 
results could be obtained. Hence, we propose the following order parameter 

(•2-7T 



1 2 
5(0, w p (</>)) 



(5) 



for the center phase transition. The novel condensate £q,^ is periodic in 
the boundary conditions and transforms like the conventional Polyakov loop 
under center transformations [10]. In contrast to the quark condensate 
{tpijj)^ the new condensate is also finite away from the chiral limit. 

Before we substantiate our proposals with numerical results we proceed 
with the investigation of the center phase transition in scalar QCD. 

Fundamentally Charged Scalar QCD 

Replacing the quarks by fundamental charged scalar fields provides a 
simpler model system due to the absence of Dirac structure. Neglecting the 
four-particle interactions the DSE for the scalar propagator read^] 

D- S \p) = Z 3 (p 2 +m 2 )-Z 3 C F g 2 TY^ J j^{p+kYD s {k)T» s {k, P] q)D^ \q) 

(6) 

with Dsip, cop) = Z$(p, uj p ) I (p^+ujp). For the scalar-gluon vertex we employ 
the model 

Ts(k, P ; q) = z, Dsl{p2 y_ D / {k2) (p + uy 

, \ A 2 Q 2 f/3oa(/^)m[g 2 /A 2 + l]y^ (7 ' 



A 2 + q 2 q 2 + A 2 \ 47T 



A related order parameter is the dual quark mass parameter M as proposed in [8]. 
For details on the derivation and the truncation scheme see e.g. |13| and references 
therein. We have omitted renormalization constants which cancel due to our scalar- 
gluon vertex model. Additionally, we apply a MOM scheme to eliminate Z3. 
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Here, the conventions for the momenta as well as the UV vertex behaviour 
is equivalent to the one in Eq.Q. In addition, we introduced a scalar 
propagator dependent function which is motivated from Ward identities in 
scalar electrodynamics and represents a generalized Ball-Chiu vertex |14j . 
The remaining part is again purely phenomenological where a dimensionless 
modeling constant d\ = 0.53 is introduced which improves the description 
of the phase transition. With the same lattice gluon propagator as in 0, 
the system can be solvecj^jand the order parameters can be extracted. 

Similar to the previous QCD case the following object serves as an order 
parameter for the center phase transition in scalar QCD 

w p (4>) 

i.e., ^s,<f> is periodic in <j) and transforms like the conventional Polyakov loop 
under center transformations [10] . 



3. Results 

In order to confirm that £q, Eq.([5]), and £5, Eq.Q, are well-defined 
order parameters we consider both QCD and scalar QCD at finite tempera- 



tures. In Fig 1(a) the ^-dependence of £q,0 and the quark condensate (ipi/j)^ 
are shown. While both quantities are periodic in [0, 2ir) and symmetric in 
(j) — >• —(f) the scalar quark dressing function B(p), evaluated at vanishing 
momenta and lowest Matsubara mode, is not periodic. Above the critical 
temperature T c = 277 MeV all three quantities melt and, in the chiral limit, 
a plateau is formed. Below T c the quark condensate is nearly constant while 
£q,0 is slightly enhanced around (f> = tt. This yields a non-zero value of £q 
below T c as compared to the dual condensate Si which almost vanishes in 



this region as shown in Fig. 1(b). Note that the temperature behaviour of 
£q below the phase transition can be further tuned by varying the vertex 
model parameters. In contrast, the order parameter £5, as defined in [6], 
considerably deviates from zero already below the phase transition, whereas 
Si and Sq stay close to zero. In order to see the expected behaviour of the 
order parameters the incorporation of all Matsubara modes is mandatory. 
For scalar QCD we obtain similar results as is demonstrated in Fig 2(a)| 



where the ^-dependence of the condensate £5^ is shown. Below T c the 
condensate is ^-independent and therefore £5 vanishes while above T c it is 
suppressed around <f> = tt and hence £5 is non- vanishing. This is displayed 



in Fig 2(b) where a distinct phase transition is observed. However, the 



results strongly dependent on the used scalar-gluon vertex model. 



For this the program CrasyDSE |15] has been used. 
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Fig. 1. Left panel: The quark dressing function £>(0, u>o(<f>)), the quark condensates 
and £<2,0, as defined in Eq.([5]), as a function of the boundary angle for 
different temperatures in the chiral limit (dashed lines: T — 273 MeV, solid lines: 
T = 283 MeV). Right panel: The order parameters Eq, Ei and E^ as defined in 
[5] as a function of the temperature. 
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Fig. 2. Left panel: Eg^ as a function of the boundary angle for three different 
temperatures and m = 1.5 GeV. Right panel: The order parameter Eg as a function 
of the temperature for three different vertex parameter values. 



4. Conclusions 



Within the Dyson- Schwinger formalism we investigated the center phase 
transition of QCD and fundamentally charged scalar QCD. A novel order 
parameter was introduced which exhibits an improved behaviour below T c 
compared to the order parameter E# proposed in [6]. Motivated by the 
strong model dependence on the used vertex further studies of the corre- 
sponding scalar/quark-gluon vertex are on-going, see e.g. 
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